home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / HQR.z / HQR
Text File  |  1996-03-14  |  4KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____HHHHQQQQRRRR((((3333FFFF))))                                                              ____HHHHQQQQRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      HQR, SHQR  -  EISPACK routine.  This subroutine finds the eigenvalues of
  10.      a REAL UPPER Hessenberg matrix by the QR method.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  hhhhqqqqrrrr((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, hhhh,,,, wwwwrrrr,,,, wwwwiiii,,,, iiiieeeerrrrrrrr))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn hhhh((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee sssshhhhqqqqrrrr((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, hhhh,,,, wwwwrrrr,,,, wwwwiiii,,,, iiiieeeerrrrrrrr))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  20.           rrrreeeeaaaallll             hhhh((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      On INPUT
  26.  
  27.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  28.      as declared in the calling program dimension statement.
  29.  
  30.      NNNN is the order of the matrix.
  31.  
  32.      LLLLOOOOWWWW and IGH are integers determined by the balancing subroutine  BALANC.
  33.      If  BALANC  has not been used, set LOW=1, IGH=N.
  34.  
  35.      HHHH contains the upper Hessenberg matrix.  Information about the
  36.      transformations used in the reduction to Hessenberg form by  ELMHES  or
  37.      ORTHES, if performed, is stored in the remaining triangle under the
  38.      Hessenberg matrix.  On OUTPUT
  39.  
  40.      HHHH has been destroyed.  Therefore, it must be saved before calling  HQR
  41.      if subsequent calculation and back transformation of eigenvectors is to
  42.      be performed.
  43.  
  44.      WWWWRRRR and WI contain the real and imaginary parts, respectively, of the
  45.      eigenvalues.  The eigenvalues are unordered except that complex conjugate
  46.      pairs of values appear consecutively with the eigenvalue having the
  47.      positive imaginary part first.  If an error exit is made, the eigenvalues
  48.      should be correct for indices IERR+1,...,N.
  49.  
  50.      IIIIEEEERRRRRRRR is set to Zero       for normal return, J          if the J-th
  51.      eigenvalue has not been
  52.         determined after a total of 30*N iterations.  Questions and comments
  53.      should be directed to B. S. Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE
  54.      NATIONAL LABORATORY
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.